Quantum 2D Heisenberg antiferromagnet: 
bridging the gap between field-theoretical and semiclassical approaches 
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The field-theoretical result for the low-T behaviour of the correlation length of the quantum 
Heisenberg antiferromagnet on the square lattice was recently improved by Hasenfratz [Eur. Phys. J. 
B 13, 11 (2000)], who corrected for cutoff effects. We show that starting from his expression, and 
exploiting our knowledge of the classical thermodynamics of the model, it is possible to take into 
account non-linear effects which are responsible for the main features of the correlation length at in- 
termediate temperature. Moreover, we find that cutoff effects lead to the appearance of an effective 
exchange integral depending on the very same renormalization coefficients derived in the frame- 
work of the semiclassical pure-quantum self-consistent harmonic approximation: The gap between 
quantum field-theoretical and semiclassical results is here eventually bridged. 

PACS numbers: 05.30.-d, 05.50.+q, 75.10.-b, 75.10.Jm 



In the last two decades, thermodynamic properties of 
the quantum Heisenberg antiferromagnet on the square- 
lattice (QHAF) have been determined by a number of 
substantially different methods. Theoretical predictions 
were compared with experimental data for real com- 
pounds [jj], as well as with numerical simulations ob- 
tained by different methods 0, 0], and also with high- 
T series expansions Despite this effort, however, a 
comprehensive picture of the subject has not yet been 
formulated. 

Much of the analysis and debate on the QHAF has 
focused on the temperature- and spin-dependence of the 
staggered correlation length ^ (T, S) . Goal of this paper 
is to show that ^ (T, S) can be expressed as 
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where the effective exchange integral Jcs{T,S) embod- 
ies quantum effects and is defined in terms of purely 
quantum spin- fluctuations, while ^c\{T/Jci) is the cor- 
relation length of the classical HAP. In particular, we 
find that Eq. , besides being the outcome of the semi- 
classical pure-quantum self- consistent harmonic approx- 
imation (PQSCHA) [3, remarkably holds also for the 
quantum field-theoretical prediction as recently im- 
proved by Hasenfratz 0. On the other hand, we show 
that the PQSCHA, when properly designed to such pur- 
pose, allows for a correct description of the low-T regime, 
via the appearance of the very same Jeff (T, S) implicitly 
entering Hasenfratz's expression. 

It is a definite suprise that a typical semiclassical ex- 
pression such as Eq. come out of a field-theoretical 
result, especially in the case of the QHAF, where semi- 
classical and field-theoretical approaches seemed destined 
to describe different regions of the (T, S) plane, with no 



possible overlap. 

The QHAF is defined by the spin Hamiltonian 



(2) 



where J>0, i—{ii,i2), d=(±l,±l), and lengths are in 
lattice units; the spin operators Si — {S^, S^, Sf) satisfy 

The correlation length ^ (T, S) is defined via the asymp- 
totic behaviour, lim^j.^^ ^G{r) cx e^l'''/^, of the stag- 
gered correlation function G{r)={—lY^^^'^ {SlSl^^.) , 
with r=(ri,r2) any vector on the square lattice. 

The classical counterpart of the QHAF corresponds to 
the hmit 5^00 with JS'^^Jd, that gives the classical 
Hamiltonian 7i=( Jci/2) Sj-Sj+d , where Sj are clas- 
sical unit vectors and Jd is the only energy scale of the 
model. 

Let us consider the right-hand side of Eq. The 
3-loop analytical expression for the classical correlation 
length, as from the field-theoretical approach Q, is 
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This expression is quantitatively correct up to T / Jd ^ 
0.4 . Above this temperature it smoothly connects to 
the exact Monte Carlo data for the classical model , so 
that the exact classical correlation length S,c\{T/ Jc\) gets 
available for all temperatures. 

As for the left-hand side of Eq. (^) , the most celebrated 
field-theoretical result is the 3-loop expression derived 
from the correlation length of the quantum nonlinear cr- 
model (QNLo-M) tSi), 
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The mapping to the QHAF consists of identifying c and 
yOg with the zero-T renormalized spin velocity and spin 
stiffness of the system 
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with Zc{S) and Zp{S) quantum renormalization coeffi- 
cients; for c and we will hereafter use the most accu- 
rate available values |lO| . as from their expansion up to 
0(S'-2) and O^S'^), respectively. 

One of the essential features of Eq. Q is the temper- 
ature-independent pre-exponential factor (e/8)(c/27rpg), 
which contrasts with the 0{T) prefactor of the classical 
Eq. lO. The 0(1) prefactor gives a purely exponential 
asymptotic behaviour ^ ~ e'^/"^, with A ^ ^-kJS"^, which 
is also the outcome of other theoretical methods such 
as the Schwinger-boson approach |lll | and the modified 
spin- wave theory [l^l- Thus the behaviour ^ ^ e^^^ has 
been generally regarded as a signature of the quantum 
character of the model. Early numerical studies sug- 
gested this signature is observed in the experimentally 
accessible temperature region of 5 = 1/2 antiferromag- 
nets |l3l| . In the last few years, however, this assumption 
has proved misleading in subtle ways, as severe difficul- 
ties arose when approaching the S> \ case by the effective 
field-theory. Finally, in Ref. it was clearly shown that 
Eq. Q holds only for temperatures low enough to ensure 
an extremely large correlation length, e.g., ^ ^ 10^ for 
5 = 1, ^ > 10^^ for S — 3/2, and generally cosmological 
correlation lengths for S > 3/2. 

In Ref . (3 Hasenfratz explained why cutoff effects, which 
are so devious for 5" = 1/2, significantly modify the cor- 
relation length for S > 1. By exploiting a direct map- 
ping between the QHAF and the QNLctM, Hasenfratz 
obtained the cutoff-corrected field-theoretical result 
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where C(T, S), defined in Eq. (14) of Ref. 7, is an integral 
of familiar spin-wave quantities over the first Brillouin 
zone. With this correction, which is the leading order 
in the spin- wave expansion for the cutoff correction, it is 
possible to obtain numerically accurate results down to 
^ > 10"' for all S. Eq. © is therefore the best available 
prediction of the field-theoretical approach. 

Our first step is to note that the function C{T, S) may 
be written as 



27rJ6'^Ci(0,-S') 
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where SC^iT, S) = ^T, S) - Cz(0, S) = 0(T/52), and 
UT,S)^1+±^^Y.^^^^(^^S)^ (8) 
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with sums over wavevectors k = {ki, in the first Bril- 
louin zone, and 7fc = (cosfci -|-cosfc2)/2 . 
From the above formulas we find 
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notice that ai is an exact coefficient, while 012 is given 
at the first order in 1/S as from Ref. 7: this is a signa- 
ture of the different approaches used in determining ^31 
and C(T, S). As for as, it contains the 3-loop correction 
terms of the field-theoretical results. 

In order to single out the relevant temperature scale in 
the exponential factor of Eq. (|lll) . one can embody the 
0{S^^) terms of the leading exponential in the argument 
of Cci(3i) • to this end we explicitly write Zp and Z^. from 
their spin-wave theory expression 10] 



Z, = Ci{0,S) + AZ, 



Ci(0,5)Co(0,^) + AZ 



where AZc and AZp are 0{S ^), and find 
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Jcff(H)(T, S) = JS^ Co(T, S) Ci(T, S) , 
a(T,5)- [(l + 5Ci+<5Co ' 
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At this level, the classical correlation length enters 
Eg. lfT^ by its 3-loop NLaM expression, Eq.(|3I). However, 
inspired by the PQSCHA result of Ref. |5|, we are led to 
generalize Eq. 1)1 3|l by replacing ^ci(3i) with the exact ^ci, 
which is numerically available at all temperatures. This 
is a fundamental step in order to get Eq. (|13|l to repro- 
duce the experimental and quantum Monte Carlo (QMC) 
data in the intermediate temperature range, a goal that 
Eq. jnj does not accomplish yet. 

In Fig.n]we show the ratio ^/^3i as a function of T/JS'^ 
for S = 5/2. The solid line is obtained by replacing ^ci(3i) 
in Ea. (|13|) with the exact ^ci- Compared with Hasen- 
fratz's expression (dash-dotted line), the agreement with 
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t = T/JS2 

FIG. 1: Ratio ^/^ai vs T/JS^ for 5" = 5/2. Dash-dotted 
line: Hasenfratz's result dotted and solid lines: Eqs. Ulljl 
and 11311 with Cci(3i) replaced by the exact ^ci; thin dashed 
line: Eq. with Jeff{H); dashed line: PQSCHA result, i.e., 
Eqs. (0 and 11611 : thin solid and dash-dotted lines: ^d/Csi s-nd 
Cci{3i)/53i, respectively (see text); the symbols are our QMC 
data. 



our QMC data greatly improves, just because the classi- 
cal correlation length is accounted for exactly. 

It is worthwhile mentioning that one could use the ex- 
act ^ci already in Eq. Hll() (dotted line in Fig. ^: such 
curve shows that the pronounced minimum in the ra- 
tio a feature that field-theoretical results do not 
reproduce, is due to substantially classical nonlinear ef- 
fects which need a non-perturbative (or exact) treat- 
ment; this is confirmed by comparing the classical ex- 
act ^ci (thin solid) and 3-loop ^ci(3i) (thin dash-dotted) 
lines, reported in Fig.^by fixing Jd to the classical limit 
{Ck—0) of Jcff(H)- The correct position of the minimum, 
on the other hand, is found by properly singling out the 
dominant temperature scale for the quantum correlation 
length, i.e., by determining the appropriate effective ex- 
change integral. To this respect, by comparing Eq. I14|) 
with Eqs. © and ((T^ . one can see that JcS{k){T, S) 
tends to for T— >0 and can be actually interpreted as 
a temperature- dependent spin stiffness Ps{T). 

Finally, we drop the 3-loop correction terms, set the 
residual a(T,S) to unity, and eventually obtain Eq. 
with the effective exchange integral Eq. H14() . The corre- 
sponding curve is shown in Fig. 1 by the thin dashed line. 
The anomalous behaviour seen as T ^ has no physi- 
cal meaning: in fact, when considering ^ as in Eq. 
with Jeff(H)j the ratio ^/Csi contains a vanishing factor 
exp(— 27r JS'^AZp/T), due to the fact that in ^31 is 
taken up to 0{S~^), while Joff(H) is only taken up to 
0(S~'^). To this respect, it is worthwhile noticing that 
Hasenfratz's expression (0 is spurious as far as the order 
in 1/S is concerned, because C{T,S) and ^31 are accu- 
rate up to 0{S~^) and 0{S~^), respectively. This origi- 



nates few other inconsistencies in the 1 /S approximation, 
essentially contained in the 0{S~'^) terms of Eq. H15|l . 
which result in slight discrepancies from the expected 
values. Nevertheless, the thin dashed line reproduces the 
behaviour drawn by QMC data in the whole temperature 
range and bridges the low- and intermediate-T regime, a 
success never scored before. 

Let us now comment on the above result. First of all, 
we notice that Jcff(H) depends just on the pure-quantum 
part of the spin-fluctuations, i.e., the difference between 
the quantum and the classical spin-fluctuations Q| (the 
signature of this being the Langevin functions Eq. ©). 
In particular, the renormalization factors Co(T^S) and 
Ci (r, S) represent the effect of the pure-quantum fluc- 
tuations of each spin with respect to its local alignment 
axis and to its nearest-neighbors, respectively. 

By considering Eq. (|14() with the low-T expansions 



SCoiT, S) 



T 



■In 



T^o 27rJS'2Ci(0,S') T 

one sees that it is the pure-quantum on-site renormal- 
ization parameter Co (T, S) that cancels the classical-like 
pre-exponential factor. This suggests that the asymp- 
totic behaviour of Eq. sets in when pure-quantum 
fluctuations of one spin relative to its nearest neighbors 
become negligible, and spins are mainly affected by on- 
site fluctuations of pure-quantum origin. 

At this point, we recall that the pure-quantum renor- 
malization coefficients Co{T,S) and (i(T,S) are charac- 
teristic of the PQSCHA, by which the correlation length 
of the QHAF in the form of Eq. ^ was first obtained 
with 



(16) 



The resemblance between Eqs. (114(1 and (|lfi|l is remark- 
able, given the fact that the PQSCHA is a completely dif- 
ferent method from that used by Hasenfratz. However, 
the PQSCHA fails to describe the regime of very low-T, 
e.g., for S=l/2 it holds only for T>0.2J [IIIEI- 

In particular, the PQSCHA has been criticized because 
the expression Eq. H16|l does not lead to the correct low- 
T asymptotic behaviour of ^. In fact, we here show that 
the PQSCHA may be properly designed to describe the 
low-T regime, thus reproducing Eq. though with c 
and only given at the first order in 1/5. This is an 
essential issue, as it means that the QHAF, even in the 
low-T limit, is a system with separable classical and pure- 
quantum aspects. 

One of the main steps in deriving the PQSCHA for 
a magnetic system is the construction of an effective 
classical spin Hamiltonian from that written in terms of 
canonical variables For the sake of clarity, we use 

(Pi Q)={Pij Qi} a-nd s={si} to indicate the set of classical 
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canonical variables and unit vectors, respectively, needed 
to describe the classical system (see Ref. for detailed 
definitions). For the partition function of the QHAF one 
finds Z= J dpdqexp{—?i{p,q)cg/T) with 

nMp,q) = o'H{§4)^e^H{s) (17) 

and d^={l + l/25)~^Ci; the further scaling p ^ 9p and 
q ^ Oq reproduces the classical-like Heisenberg Hamilto- 
nian H{s)=J{S+l/2)'^9'^ J^id ^i-^i+d , thus defining Joff 
as in Eq. (|f 6|l . When the correlation functions are consid- 
ered, different scaling laws come into play, viathe integral 

G(r) cx Jdp dq g^s{r; p, q) e-^-'t^P'-?)/?^ (^^g) 

where (/cff scales according to 

5eff(r-;p, q) = et g{r; ^, ^) 9^. g{r; s) (19) 

with g{r; s)=N^^J2i^i'^i+r', &i. is a pure-quantum 
renormalization coefficient defined in Ref. 0. As both 
'H{p,q) and g{r;p,q) are biquadratic functions, the dif- 
ferent scalings (with or 9r) conflict. At intermediate T, 
thermodynamics is governed by many different spin con- 
figurations, and the effective Hamiltonian appearing in 
the Boltzmann factor in Eq. (|18|l dominates (i.e., scaling 
with 9 is preferable 0): this gives Eq. H16|l . At low- 
T the quartic terms in the effective Hamiltonian do not 
matter, and the scaling conflict is settled: therefore, the 
best evaluation of the integral ((TH|l is made by scaling 
with 9r, so that the prefactor 6*^ in Eq. (|17|l is replaced 
by 9^9^. As for increasing \r\ the coefficient 9r rapidly 
converges to {l+l/2S)-^Co: the PQSCHA effective ex- 
change for the low-T correlation length reads JS^ Co Ci : 
consistent with Eq. (|14|) . 

A quantitative comparison between Hasenfratz's and 
PQSCHA results should also take into account some im- 
portant details concerning the spin- wave dispersion rela- 
tion used to evaluate the pure-quantum renormalization 
factors Eqs. (|SJ). We just mention that in order to avoid 
ambiguities due to the ordering prescription one should 
definitely prefer the zero-T renormalized frequencies H10() 
to the bare frequencies originally used by Hasenfratz. In 
fact, the PQSCHA method 0, Q uses even more re- 
fined temperature- and configuration-dependent renor- 
malized frequencies; treating them within the so-called 
low- coupling approximation yields a more accurate finite- 
T expression, which is crucial for quantitative accuracy 
in evaluating thermodynamic quantities . 

In this paper we have shown that, when cutoff effects 
are properly taken into account, the field theoretical ex- 
pression for the correlation length of the QHAF analit- 
ically contains its classical counterpart, according to a 
simple equation that gives the former in terms of the lat- 
ter with an effective exchange integral. Such equation 
surprisingly holds at all temperatures. Once the exact 



classical correlation length is available, the above result 
is used to make clear that the main features of the quan- 
tum correlation length at intermediate temperatures are 
due to essentially classical non-linear effects, which can- 
not be taken into account by perturbative approaches. 
Moreover, the effective exchange integral is seen to de- 
pend on the same pure-quantum renormalization coeffi- 
cients defined by the PQSCHA, according to an expres- 
sion which is very similar (equal) to that found by the 
latter approach in its standard (low-T) version. 

The idea that an effective exchange integral may be 
uniquely defined for all thermodynamic quantities of the 
QHAF, as predicted by the PQSCHA, is suggestive and 
deserves, in our opinion, further investigation by the 
field-theoretical community. 
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